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We study the dynamics of a lattice of coupled nonidentical Fitz Hugh-Nagumo system subject to 
independent external noise. It is shown that these stochastic oscillators can lead to global synchro- 
nization behavior without an external signal. With the increase of the noise intensity, the system 
exhibits coherence resonance behavior. Coupling can enhance greatly the noise-induced coherence 
in the system. 
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The study of coupled oscillators is one of the funda- 
mental problems with applications in various fields |Q. 
Mutual synchronization of the oscillators is of great in- 
terest and importance among the collective dynamics of 
the coupled oscillators. The notation of synchronization 
has been extended to include a variety of phenomena in 
the context of interating chaotic oscillators, such as com- 
plete synchronization H, generalized synchronization [ 3 1 , 
phase synchronization and lag synchronization [3|. 
Also, synchronization phenomenon has been studied in 
stochastic systems. Stochastic resonance can be under- 
stood from the view point of frequency locking and phase 
synchronization of the noise-induced motion to the ex- 
ternal signal M. Due to this synchronization, stochas- 
tic bistable elements subject to the same periodic signal, 
when appropriately coupled, can display global synchro- 
nization to the external periodic signal. This synchro- 
nization has the effect of enhancing the stochastic reso- 
nance in the array || . 

Noise-induced coherent motion has been a topic of 
great interest recently. For a system at a Hopf bifur- 
cation point, an optimal amount of noise can induce the 
most coherent motion in the system ||. More generally, 
noise-induced coherent motion has been demonstrated in 
a variety of excitable systems and system with 

delay . Although described by different terms such as 
stochastic resonance without external periodic force, au- 
tonomous stochastic resonance or coherence resonance, a 
common feature of this type of systems is the increased 
coherence of the motion and resonant like behavior of the 
coherence induced purely by noise without an external 
signal. An interesting question about this type of systems 
is how these stochastic elements behave when coupled. 
Do the elements subjected to independent noise display 
synchronization? and can the synchronization enhance 
the coherence in the system motion? A recent study 
demonstrated noise-induced global oscillation and reso- 
nant behavior in a subexcitable media Jl4| ]. In Ref. Jl5| , 
it is shown that two interacting coherence oscillators can 
be synchronized. However, whether the coherence of the 
motion can been enhanced by the interaction between 



the elements are not clear. 

This paper studies these problems by the following N 
coupled nonidentical Fitz Hugh-Nagumo (FHN) neurons, 
a simple but representative model of excitable system and 
nerve pulses jill, 



Xi- — -Vi + g\Xi+i + Xi-i 



y i = Xi + a. L + D£i(t), 
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where Oj is a parameter of the ith element. For a single 
FHN model, if \a\ > 1, the system has only a stable fixed 
point, while \a\ < 1 a limit cycle. The system with fixed 
point dynamics (|a| slightly larger than one) is excitable 
because it will return to the fixed point only after a large 
excursion ( "near limit cycle" ) when perturbed away from 
the fixed point. To make the study more general, we sup- 
pose di is not the same for the elements, but has a uni- 
form distribution a, G (1,1.1). This implementation of 
the model is physically significant because many physical 
systems are diffusively coupled, and nonidentity is more 
natural in physical situations. With nonidentical aj, the 
uncoupled elements will have different response to the 
same level of noise £j(i), e.g., different average firing fre- 
quency. The Gaussian noise is uncorrelated in different 
elements, i.e., (£,i(t)£,j(t')) = SijS(t — t'). Periodic bound- 
ary condition is employed in our studies. The parameter 
of the time scale is fixed at e = 0.01. 

To characterize synchronization behavior in the lattice 
of nonidentical excitable systems, we introduce the phase 
of the elements Iml 



t - Tfc 



27rfc, 



Tk 
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where Tk is the time of the fcth firing of the element 
defined in simulations by threshold crossing of Xi(t) at 
x = 1.0. The quantity 



(4) 



measures the phase synchronization effect of neighbour- 
ing elements H. A spatiotemporal average of Si, i.e., 
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gives a measure of the degree of phase synchronization 
in the coupled system. For completely unsynchronized 
motion S ~ 0.5, while for globally synchronized system 

To measure the temporal coherence of the noise- 
induced motion, we examine the distribution of the pulse 
duration Tk = Tk+\ — Tk- For a single element subject to 
noise, the distribution P(T) has a peak at a certain value 
of Tk and an exponential tail at large values |lC| . A mea- 
sure of the sharpness of the distribution, for example, 

R= (T,)/VVar(T fc ), (6) 

which can be viewed as signal-to-noise ratio, provides an 
indication of the coherence of the firing event. Biologi- 
cally, this quantity is of importance because it is related 
to the timing precision of the information processing in 
neural systems [jl6). For a single element, it has been 
shown that R possesses an optimal value at a certain 
level of noise M] . Here we adopt the same measure of co- 
herence, with the distribution P(T) constructed by pulse 
duration of all the N elements during a long enough pe- 
riod of time. 

Typical behaviors observed in numerical simulations of 
the system are in the following: 

(1) Without external noise, each element comes to a 
fixed point state and no firing takes place. The fixed 
points are slightly different due to the nonidentity. 

(2) When strong enough noise is added to the origi- 
nally quiescent system, the system is excited. The firing 
process of the elements may be in phase synchronization 
if the coupling is strong enough. 

(3) The coherence of the temporal and spatial pattern 
of the firing process is greatly enhanced in a region of 
noise and coupling strength. 

The main features of the system are illustrated in Fig. 
1, where the results of R and S in the parameter space 
(log 10 (<?), log 10 (.D)) are shown for a lattice of N = 100 el- 
ements. For a fixed coupling g, R increases first with the 
noise level D, reaches a maximum, then decreases again, 
showing the typical resonant behavior without an exter- 
nal signal. In general, for a stronger coupling, a higher 
level of noise is needed to excite the system. Similarly, 
for a fixed noise level D, R increases with increasing g 
until it reaches an optimal value; after that, it decrease 
again. 

From Figure 1, one can observe several dynamical 
regimes in the systems. For very weak coupling (g < 
10~ 2 ), the firing of the elements are essentially indepen- 
dent, because a noise-induced firing of an element cannot 
excite its neighbouring lattices. Due to the independent 
firing, the phase difference has a uniform random distri- 
bution on (0, 27r), resulting S ~ 0.5. In this region, with 
the increase of D, the temporal behavior of the system 
display coherence resonance similar to that observed in 



a single element. The enhancement of the coherence by 
the noise is not very pronounced. A maximal R rs 4.5 
is found around D « 10 _1 . A typical spatiotemporal 
pattern of x% and Si for weakly coupled elements subject 
to relatively weak noise is shown in Fig. 2(a). Both the 
spatial and temporal behaviors are quite irregular. Each 
element has its individual average firing frequency due to 
the nonidentity (smaller frequency for larger aj). 

With the increase of the coupling strength, the system 
becomes sensitive to weak noise because the firing events 
induced by noise now become the source of excitation 
of the neighbouring elements. This mutual excitation 
enhances the coherence of the motion in the coupled sys- 
tem, as indicated by increasing R and decreasing S. The 
lattice displays clusters of synchronization. The clusters 
breaks and reunites during the evolution, so that each 
element has slightly different firing frequency. Typical 
behavior of this partial synchronization regime is shown 
in Fig. 2(b). However, if the external noise is strong 
enough so that the noise dominates over the coupling, 
firing of each element is governed mainly by its individ- 
ual noise, and synchronization clusters cannot survive, 
as seen in Fig. 2(c). The coherence of the temporal be- 
havior is relatively low because quite large noise deforms 
greatly the "near limit" cycle. Note that firing frequency 
is not locked, but the difference is not as pronounced 
as in weak noise region (Fig. 2(a)). The next regime 
where R takes large values (R ~ 18) while S is very 
close to zero, is the most interesting, because the system 
performs quite regular motion globally, as seen in Fig. 
2(d). All elements are locked to a relatively large firing 
frequency, and the distribution of the pulse duration be- 
comes very sharp. After that, with stronger coupling, the 
system keeps global synchronization, however, the tem- 
poral behavior becomes irregular again, as indicated by 
decreasing R. This can be understood qualitatively from 
the global dynamics X = (xi)^ and Y — (Yi)ry) with () jy 
denoting average over the lattice. From Eqs. (1-2), we 
get approximately 

X 3 

(X = X r 2 X - y, (7) 

Y = X + a Q + -^Lm, (8) 

where r 2 = ((xi — A) 2 )^ is the fluctuation level of lo- 
cal dynamics and higher order terms of this fluctuation 
are ignored; ao = (ai) n- The summation of the indepen- 
dent noise is still a Gaussian noise £(t), but with a weaker 
strength D / y/N. For strong enough coupling, the system 
achieves global synchronization, Xi « X, and r 2 rs 0, and 
the coupled system can be viewed as a single element sub- 
ject to a white noise with strength D/y/N. For a single 
element, a rather weak noise may not excite the system, 
or the excited motion is quite irreg ular 0]. This ex- 
plains the globally synchronized but irregular motion of 
the lattice, as shown in Fig. 2(e). 
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The locations of the above five representative dynam- 
ical regimes shown in Fig. 2 have been indicated by the 
black dots in Fig. 1. These regimes are typical for dif- 
ferent size of the lattice. For a larger lattice, the regime 
of global regular motion (R large and S ~ 0) is wider in 
the parameter space. 

Now let us discuss the phenomenon of array enhanced 
coherence resonance. For a fixed size of lattice and a 
certain value of coupling g 7 there is an optimal level of 
noise at which the system has a maximal value R m ax- 
Figure 3 shows R m ax a s a function of g for different size 
of the lattice. The dashed line represents R ma x for a 
single element with a = 1.05. As seen from this figure, 
in the weak coupling region, R max of the coupled lat- 
tice is lower than that of a single element, because the 
nonidentity makes the distribution of the pulse duration 
broader. However, in some intermediate coupling region, 
even only two coupled elements can enhance the coher- 
ence of the temporal behavior, even though they are non- 
identical. The enhancement is larger for larger N. For 
large enough lattice, R m ax seems to be saturated, but 
the region of coupling in which R max takes large values 
is broader for larger lattice. As pointed above, for strong 
enough coupling, the coupled system tends to act as a 
single element, and Rmax converges to that of the single 
element at large value of g. Clearly, for a larger lattice, 
stronger coupling is needed to make the whole lattice act 
as a single element due to the local diffusive coupling, 
resulting in a slower convergence of R ma x to that of the 
single element. 

The properties observed here have some similarity to 
that in coupled bistable stochastic resonance oscillators 
subject to the same periodic signal The difference 
is that in Refs. || and ]l4|j , global synchronization is 
observed only at an optimal level of noise, while in the 
present system, global synchronization occurs for strong 
enough coupling if the noise level is not too high. The 
global motion can be regular or irregular. Here we should 
emphasize that there is no periodic signal in our system, 
and the coherent motion is purely induced by noise and 
enhanced by coupling. 

To conclude, we demonstrate global synchronization 
and array enhanced coherence in a lattice of locally cou- 
pled, nonidentical FHN model neurons. The results are 
similar for identical neurons. Interaction between the 
elements not only renders synchronization of the firing 
process induced by independent noise, but also improves 
the coherence of the noise-induced motion greatly. The 
phenomena demonstrated may be of importance in neu- 
rophysiology, where background noise may play an active 
role to increase the order and timing precision of a large 
ensemble of interacting neurons in biological information 
processing. 
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FIG. 1. Signal-to-noise ratio R (a) and degree 

of phase synchronization S (b) in the parameter space 
(log 10 (i?), log 10 (D)) of a coupled lattice of nonidentical FHN 
neurons with N = 100. 

FIG. 2. Five typical dynamical regimes in the parameter 
space. The upper panels show the spatiotemporal structure 
of Xi and Si. The time step is 0.2. The lower panels show the 
average firing frequency of the lattice and the distribution of 
the pulse duration. The location of these representative dy- 
namical behaviors are shown by black dots in the parameters 
space in Fig. 1. The parameters are (a) g — 0.005, D = 0.02; 
(b) g = 0.03, D = 0.025; (c) g = 0.02, D = 0.25; (d) 
g = 0.25, D = 0.07; and (e) g = 1.0, D = 0.04. 
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FIG. 3. Illustration of array enhanced coherence reso- 
nance. The maximal value of R is plotted as a function of g 
for various sizes of the lattice. 
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